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ABSTRACT 

Aims. We study selected properties of solar energetic particle (SEP) events as inferred from their associated radio emissions. 
Methods. We used a catalogue of 115 SEP events, which consists of entries of proton intensity enhancements at one AU, with 
complete coverage over solar cycle 23 based on high-energy (~ 68 MeV) protons from SOHO/ERNE. We also calculated the proton 
release time at the Sun using velocity dispersion analysis (VDA). After an initial rejection of cases with unrealistic VDA path lengths, 
we assembled composite radio spectra for the remaining events using data from ground-based and space-home radio spectrographs. 
We registered the associated radio emissions for every event, and we divided the events in groups according to their associated radio 
emissions. In cases of type Ill-associated events, we extended our study to the timings between the type III radio emission, the proton 
release, and the electron release as inferred from VDA based on Wind/3DP 20-646 keV data. 

Results. The proton release was found to be most often accompanied by both type III and II radio bursts, but a good association 
percentage was also registered in cases accompanied by type Ills only. The worst association was found for the cases only associated 
with type II. In the type Ill-associated cases, we usually found systematic delays of both the proton and electron release times as 
inferred by the particles’ VDAs, with respect to the start of the associated type III burst. The comparison of the proton and electron 
release times revealed that, in more than half of the cases, the protons and electrons were simultaneously released within the statistical 
uncertainty of our analysis. For the cases with type II radio association, we found that the distribution of the proton release heights 
had a maximum at ~2.5 Rq. Most (69%) of the flares associated with our SEP events were located in the western hemisphere, with a 
peak within the well-connected region of 50° -60° western longitude. 

Conclusions. Both flare- and shock-related particle release processes are observed in major proton events at >50 MeV. Typically, the 
protons are released after the start of the associated type III bursts and simultaneously or before the release of energetic electrons. Our 
study indicates that a clear-cut distinction between flare-related and CME-related SEP events is difficult to establish. 
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1. Introduction 


Solar energetic particle (SEP) events , are among the most impor¬ 


tant e l ements of space weather (e.g. Reame^ 1999t IVainio et al.l 

It 


2009t l^palswamv et al.ll2008l : IValtonen 2011 ). The energy of 
SEPs usually range from ~10 keV/nucleon up to some GeV, 
the events last from several hours to a few days, and their fluxes 
can rise above the background by many orders of magnitude. 
They accompa ny solar flares and coronal ma ss ejections (CMEs) 
(cf. reviews by lReameslll999l: ICliverllToOOl) . and both types of 
activity provide possibilities for SEP-related particle accelera¬ 
tion through magnetic reconnection and large-scale coronal dis¬ 
turbances such as shock waves. An important issue is whether 
both magnetic reconnection and shock wave acceleration can si¬ 
multaneously accelerate particles or whether one of them domi¬ 
nates the particles energization process. The scenario of diffusive 
shock acceleration seems to be a plausible source of low-energy 
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protons up to some MeV; however, it is unlikely that the shock 
mechanism alone applies to very high-energy particles from the 
Sun, in the relativistic range. 

Work in the 1990s indicated that the SEP events could be 
separated into two distinct classes based on the properties of the 
associated SXR flare, correlations with radio bursts, the pres¬ 
ence or absence of a CME, the abundances, and charge states. 
Th e SEP events were classified as ehher impulsive or gradual 
(cf. lReames|[l999l: lKallenrod3l2003t iReamesi 120131) . Impulsive 
SEP events may originate from so lar flares, usually small flares 
(C-M-class; see iNitta et al]l2006h . last less than one day, haye 
relatiyely small peak fluxes, and haye b een associated with type 
III radio bursts (iReames & Stonel[l986l) . Gradual SEP eyents are 
thought to be accelerated in CME-driyen coronal and/or inter¬ 
planetary shocks. Compared to impulsiye SEP eyents, gradual 
SEPs haye long durations (i.e. days), relatiyely high peak flu xes, 
and haye long been associated with type II radio bursts dLinl 
II970I: ISyestka & Pritzoya-Syestkoyaill974l : ICliyer et"^ l2004l) . 
Although this classification has been useful, we point out that 
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both flares and CMEs are associated with almost all SEP events, 
so it is often difficult to distinguish the particle accelerator un¬ 
ambiguously. 

Radio emission in the solar corona and interplanetary (IP) 
medium provides important information about the acceleration 
and propagation of solar energetic particles as well as shocks. 
The different types of radio bursts in the dynamic radio spec¬ 
trum, are usually interpreted a s signatures of electron beams 
(type Ills: iLin et alj l98^ 1^^ ), electrons a ccelerated in shock 


waves (type IIs, 


Nelson & Melrose! 1 19851 : ICairns et^ 120031 : 


iGanse et alj|2012|) . or a s electrons confined (type IVs) in closed 
loop structures (see e. 2 . iNindos et al.ll2008h . 

So far various studies have investigated the association be¬ 
tween t he type III, II, and IV radio bursts and SEP events in 
general. iLiril (Il970li found that SEP events from 1964 to 1967 
were mostly accompanied (70%) by type III and II radio bursts; 
subsequently [Svestka & Eritzova-Svestkoval (Il974li . using type 

II burst records and SEPs from 1966 to 1968, concluded that al¬ 

most a ll type I I radio bursts were associated with solar proton 
events. iKahleij (Il982h suggested that only about half of metric 
type II bursts with sources in the western solar hemisphere were 
linked to SEPs. Even though coronal and IP shock waves are 
well established to be the principal accelerators for the largest 
obs erved SEP events, not all c oronal shocks are "SEP effective" 
(see lGopalswamv et al.ll200^ . ICane et alJ (l2002h concluded that 
the proton events detected above ~20 MeV (Jan. 1997-May 
2001) were always accompanied by type III bursts, by CMEs 
and, for events originating not too far behind the limb, by flare 
signatures in Her and soft X-rays. Even though ICane et alJ (l2002h 
reported that all the SEP events are accompanied by CMEs, their 
study focussed only on type III radio association and a possible 
type II occurrence was not pres ented. _ _ 

The statistics presented in ICliver etTP (l2004l) . from 1996 
July through 2001 June, supported the widely held view that the 
SEPs (~20 MeV) observed near Earth were associated with coro¬ 
nal shock waves and especially with metric typ e II bursts that 
have d ecametric-hectometric (DH) counterparts. ICliver & Ling! 
(l2009ll separated their SEP sample into gradual and impulsive 
events and found that gradual events are always accompanied 
by DH-type II radio bursts, impulsive events had no associa¬ 
tion with DH-type II (except for one case), and in all cases 
an associated type III burst was observed. Also type III bursts 
with associated DH-type II bursts were ten times more likely 
to have large associated SEP events above 30 MeV (also see 
lGot)alswamy _&Ma^£ 20101) . In the most recent study pre- 
sented bv iMiteva et al.l ( 20131) . SEP events, either gradual or im¬ 
pulsive, were found to have the highest association rate with type 

III radio bursts and a lower association with type II bursts. 

The comparison of the release times, between ions, electrons, 
and electromagnetic emission can, in principle, be used to iden¬ 
tify the most relevant acceleration processes. Radio waves orig¬ 
inate from energetic electrons, thus any comparison with ions in 
space is indirect. So far, energetic electrons have been studied 
with respect to their associate d radio emissions and especially 
with t ype III radio bursts ( e.^. iHaggertv & Roeloj 200^ ICand 
I 2 QO 3 I) . [&ucker et ^ d 1 999l) and iHaggertv & Roelof (l2002h con¬ 
cluded that the release of electrons into interplanetary space of¬ 
ten starts several minutes after the first radiative signatures of 
electrons in the corona, assuming negligible interplanetary prop¬ 
agation effects. The late electron release can be ascribed to: (1) 
the action of a different acceleration process for the escaping 
electrons, i.e. the delay is due to electron acceleration at a coro¬ 
nal shock (ISimnett et alj200^: (2) transport effects during inter¬ 
planetary propagation (ICandl2003h : (3) delayed electrons origi¬ 


nating at a reconnection site behind the CME shock front at re- 
gio ns of magnetic res tructuring linked to the CME development 
(see lMaia et al.ll200^ . 

Compared with proton release time, relativistic electrons 
were found to be released with a signif icant delay in many cases 
deliver et al.lll^ . lCliver et all (Il982h concluded that 100 keV 
electrons were released first, followed by the 2 GeV protons 
within 5 minutes, and 1 MeV electrons following at least 5 min¬ 
utes later. This apparent delay in proton-electron acceleration 
may be due to several causes, such as (1) selective acceleration, 

(2) cross-field transport or perpendicular particle diffusion, and 

(3) particle trapping. 

In this article we study selected properties of high-energy 
SEP events, especially focussing on the association of the SEP 
release time, as inferred by the velocity dispersion analysis 
(VDA), with transient solar radio emissions recorded by space 
and ground-based radio spectrographs. We have also thoroughly 
examined the time difference of proton release with respect to 
the escape of keV electrons into space, derived from both the 
electron release as traced by DH type III bursts and the elec¬ 
tron release as derived from VDA. Eor the cases with type II 
radio association, we estimated the proton release heights from 
the height-time profile of the CME leading edge. 


2. Data selection and analysis 

2 . 1 . Proton event selection 

Eor the purpose of our analysis, we used a list of 115 high-energy 
solar proton events as record ed by the Energet ic and Relativistic 
Nuclei and Electron (ERNE: iTorsti et al.l 19951) instrument of the 
University of Turku, on the Solar and Heliospheric Observatory 
(SOHO) from May 1996 until the end of 2010 extending over a 
fu ll solar cycle . The c atalogue of these SEP events was presented 
in IVainio et al.l (1201 3l) and consists of entries of proton intensity 
enhancements based on the systematic scan of proton intensities 
in the energy range of 54.8-80.3 MeV. To register a proton in¬ 
tensity enhancement as a SEP event, we used the criterion that 
the 1 minute average intensity in the 54.8 - 80.3 MeV ERNE pro¬ 
ton channel should be enhanced by a factor of ~ 3 above the quiet 
time background. In Eig.[T]we present a representative case, the 
15 January 2005 SEP event. 

The energy range of 54.8-80.3 MeV in the SEP selection 
criterion was purposely chose n well above the typi cally used 
>10 MeV proton channel {e.g. iLaurenza et al.1l2009l) . available 
from GOES satellites, to identify individual events in cases in 
which large SEP events mask smaller events at low energies 
when they follow each other in quick succession. At high en¬ 
ergies these cases can be better distinguished because of a more 
rapid fall of the proton intensities. However, the ERNE High En- 
ergy Detector has a very large geometri c factor (~ 20 -40 cm^ sr; 
IValtonen et al.|[l9^ iTorsti et Sill 9951) often causing saturation 
during yery large SEP eyents. Therefore, eyents following larger 
eyents before the recoyery of the instrument may be missed. 
These cases, among others, are the eyents of 17 and 20 January 
2005 following the eyent of 15 January 2005 (see Table 1 of 
ICane et al.l2010i list of SEP eyents obseryed from 1997 to 2006, 
at energies aboye 25 MeV), which were not included in our cat¬ 
alogue. 

2.2. Estimate of proton release time 

We determined the proton onset time with the so-called Poisson- 
CUSUM method in 20 proton energy channels of SOHO/ERNE 
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SOHO/ERNE 1-min. Average: 15 / 01 / 2005 



(1.49- 1.78] 1.63M6V/n 
(1.78-2.16] 1.97 MeV/n 
(2,16-2.66] 2,41 MeV/n! 
(2.66 - 3.29] 2,98 MeV/n: 
(3.29-4.10] 3.70 MeV/n, 
(4.10-5.12] 4,61 MeV/n: 
(5.12-6.42] 5.77 MeV/n 
(6,42 - 8.06] 7.24 MeV/nl 
(8.06- 10.1] 9.09 MeV/n: 
(10.1 - 12.7] 11.4M6V/n' 
(13.8- 16.9] 15.4MeV/ni 
(16.9-22.4] 18.9M6V/n: 
(20,8 - 28,0] 23.3 MeV/n. 
(25.9 - 32.2] 29.0 MeV/n: 
(32.2 - 40.5] 36.3 MeV/n • 
(40.5 - 50.8] 45.6 MeV/nl 
(50.8 - 67.3] 54.0 MeV/n: 
(54.2 - 79.2] 67.5 MeV/n. 
(79.2- 114] 94.0 MeV/ni 
(111 - 140] 116MeV/n : 


00:00 06:00 12:00 18:00 00:00 


Velocity Dispersion Anoiysis 



20050115 Wind/Waves + ART. Dynamic Spectrum 
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Fig. 1: Proton Event of 15 January 2005: Top: proton intensity as 
recorded from SOHO/ERNE for 20 proton energy channels be¬ 
tween 1.58 MeV and 131 MeV. Middle: proton velocity disper¬ 
sion analysis from the SOHO/ERNE data. We mark the proton 
onset times in the individual energy channels with blue crosses, 
the horizontal error bar represents the uncertainty of the parti¬ 
cles’ inverse velocity, which is determined from the width of ev¬ 
ery energy channel and the vertical error bar represents the un¬ 
certainty of the onset time determination derived from the Monte 
Carlo analysis. With the straight red line we present the linear 
ht to the onset times using the maximum likelihood estimation 
approach. Bottom: composite radio spectra of Wind/Waves and 
ARTEMIS-IV. The vertical blue solid line and dashed lines de¬ 
note the proton release time derived from the VDA and its uncer¬ 
tainty, respectively; the red line denotes the onset time of 54.8- 
80.3 MeV protons. 


between 1.58 MeV and 131 MeV for each SEP event. Cumula- 
tive sum (C USUM) quality control schemes were proposed by 
iPagel (Il954li and are widely used in other fields, especially in 
medical research. A CUSUM scheme cumulates the difference 
between an observed and a reference value. If this cumulation 
equals or exceeds a threshold, then an out-of-control signal is 
given at the exact moment when the process transition to the ob¬ 
served value has occurred. 


Based on the dehnition of the CUSUM schemes, the on¬ 
set time of an event, at a given energy channel, is the instant 
when the signal has a systematic change of its mean value com¬ 
pared to the pre-event background mean. Usually, to consider 
that a transition in the signal mean value has occurred we use 
the two-sigma-shi ft criterion dHuttunen-Heikinmaa et al.1 120051 : 
IVainio et alJl2013h . Previous studies have used other methods 
for the determinatio n of the particles’ onset times: for example 
iKrucker et al.l (1 19991) used a She whart control-chart sc heme with 
an upper control limit at 4 cr and iMasson et alJ (12012h presented 
a method in which the onset times were determined by reference 
times in the background-normalized rise time prohles at differ¬ 
ent energies, which works similar to a cross-correlation analysis. 
In general, the CUSUM schemes have been shown to be very ef¬ 
ficient in detecting small shifts in the mean of a process; they are 
effective even when it is desired to detect shifts in the mean that 
are 2 cr or less. 

In each energy channel, we combined the CUSUM method 
with a novel technique based on a statistical analysis of the pa¬ 
rameters (mean and standard deviation) used by CUSUM to es¬ 
timate the uncertainties of the onset times. In the cases where the 
assumptions about the statistics of the sequence of random vari¬ 
ables are not valid, namely, the random variable does not have 
a normal probability density function and the monitored process 
is not stationary, then the selection of the time interval where the 
background mean and the standard deviation is calculated can 
produce a signihcant uncertainty in the estimated onset times. 

We estimate the uncertainty of the determined onset times 
using a Monte Carlo analysis of the CUSUM parameters, which 
are calculated inside the pre-event time interval. We separate the 
initial time interval into randomly distributed sub-domains of 
variable sizes and we calculate the signal parameters in each it¬ 
eration. Then, we estimate an onset time for every iteration using 
the CUSUM method. If the assumptions about the signal statis¬ 
tics are not valid, then different signal parameters are calculated 
in every iteration, and therefore, the onset time estimation with 
the CUSUM method could be different for every parameter set. 
In most of the cases this process produces a collection of esti¬ 
mated onset times that are normally distributed around a certain 
fixed value. The maximum of the distribution corresponds to the 
proton onset time with the highest probability and its width to 
the onset time uncertainty. 

The velocity dispersion analysis (VDA) of the time when the 
first particles arrive at the detector is a method that is frequently 
used to estimate the release time of the energetic particles and 
their interplanetary travel path length. This method relies on the 
fact that the less energetic particles arrive at Earth later than the 
more energetic particles. The VDA is based on the assumptions 
that the first arriving particles of all energies are simultaneously 
released from the same acceleration site in the solar corona and 
their propagation is scatter free in the interplanetary medium. We 
should also add as a necessary condition to the above-mentioned 
assumptions that the energetic particles have an energy density 
lower than that of the local magnetic held to be effectively con- 
hned during their propagation in magnetic Hux tubes without sig¬ 
nihcant cross-held transport. In AppendixlAlwe calculate the en¬ 
ergy density for a selected, very energetic event to prove that the 
energy density of the energetic protons is negligible compared to 
the magnetic energy density. 

The velocity dispersion equation at 1 AU can be written as 

trei = tonset{E) - s ■/?^'(£') -H 8.33 min., (1) 

where tonsetiE) is the observed onset time of the hrst arriving par¬ 
ticles at 1 AU, hei is the release time of protons with energy E 
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at the Sun, s is the apparent path length travelled by the parti¬ 
cles, and I3^\E) = ^ where u(E) is the relativistic velocity of 
protons with energy E. Erom the linear relationship between the 
arrival times of the first particles measured on Earth tonsetiE) and 
the jS '(£) the solar release time can be estimated from the in¬ 
tersection at /1(E)“' = 0 and the interplanetary distance travelled 
by the energetic particles is estimated from the slope. 

We constructed the VDA plots from the onset times that were 
determined by the CUSUM method at different energies and then 
we plotted them versus fE^iE). We also included the uncertainty 
of the onset time determination in every energy channel and the 
uncertainty of the corresponding jS '(£) value, which is deter¬ 
mined from the width of every energy cha nnel. Eor the linear 
fit, we used the method of lYork et al.l (l2004l) , which applies the 
maximum likelihood estimation approach for the determination 
of the slopes, the intercepts, and the standard errors of a number 
of data sets. This method produces a weighted least-squares fit 
of a straight line to a set of points with error in both coordinates. 

The result of the VDA for the 15 January 2005 event is pre¬ 
sented in the middle panel of Eig. [T] We mark with blue crosses 
the onset times and their uncertainty, derived from the Poisson- 
CUSUM method and the Monte Carlo method, respectively. The 
red line corresponds to the linear fit of the onset times derived 
from the maximum likelihood estimation method. Erom Eq. [T] 
we calculated the apparent path length and the proton release 
time at the Sun (the 8.33 minutes of travel time appear in Eq.[T] 
to enable the comparison between the arrival times of protons 
and the electromagnetic radiation at 1 AU). In the case of Eig. [1] 
the release time was determined to be at 06:23 UT with an un¬ 
certainty of +3.0minutes and the apparent path length 1.67 AU 
with an uncertainty of ±0.08 AU. 

Usually, the onset times are strongly affected by fluctuations 
of intensities in the individual energy channels, which may intro¬ 
duce a distortion and produce a large error both in the apparent 
path length and the release time. In our analysis, we considered 
that when the apparent path length values were outside the range 
of 1 - 3 AU, the VDA and the related release times should not be 
trusted; thus, 32 out of 115 events were rejected from our anal¬ 
ysis. Nominally, the path length should be close to the length of 
the Parker spiral from the Sun to the Earth {i.e. ~1.15 AU, for 
solar wind velocity i^s w - 440km/s). However, several velocity 
dispersion studies {e.g. iKrucker & Linll2()()(il : IVainio et al.ll2013l) 
have shown significant deviations between the path lengths that 
were obtained from the VDA and the Parker spiral length. These 
deviations can arise, in cases of larger path lengths (> 1.5 AU), 
when ener getic particles are subject ed to strong interplanetary 
scattering (iLintunen & Vaini^ |2004 path lengths from simu¬ 
lated data) or when the interplanetary magnetic field is dis¬ 
turbed (i.e when it has a significant component out of the ecliptic 
plane). Moreover, deviations fo r path lengths < 1.0 AU can arise 
from instrumental effects (e.g. iHaggertv & Roel(^l2002l 120031 : 
iKrucker et alJl200ll in the case of energetic electrons). Another 
reason is that the hypothesis used in VDA that particles at all en¬ 
ergies are simultaneously released at the same distance from the 
spacecraft may be wrong. 


Table 1: Overview of the Radiospectrographs. 


Radio 

Spectograph 

Observing 
Time (UT) 

Erequency 

(MHz) 

Temporal 

Resolution 

Wind/Waves 

All Day 

0.02-13.825 

1 min 

Artemis-IV 

06:00-15:00 

20-650^ 

1 sec 

Culgoora 

20:00-07:00 

18-1800 

3 sec 

RSTN 

All Day 

25-180 

3 sec 


Before July 2002 the frequency range was 110-650 MHz. 


olution from 20kHz to 13.825 MHz, from the two WAVES 
receivers known as RADI (20 kHz - 1.04 MHz) and RAD2 
(1.075-13.825MHz). The lower-coronal parts in our compos¬ 
ite spectra consisted of data from the followi ng ground instru¬ 
ments (radio spectrographs): A RTEMIS-IV (ICaroubalos et al.l 
1200 U l ^ntogeorgos et alJ 20061) . Culgoora solar radio spectro¬ 
graph Solar Telescope Network 

(RSTNT TGuidke et al.lll98l1) . We occasionally used Nancay De¬ 
cameter Array (DAM) quick-look spectra from the radio survey 
projecfl. We used ARTEMlS-IV data from the sweep frequency 
receiver ASG with a temporal resolution of 0.1 second, which 
we integrated over 1 second. An overview of the radio data we 
used is presented in Table [T] 

Eor the event of Eig. [1] we present the associated composite 
radio spectrum from Wind/Waves and ARTEMlS-IV on which 
we have marked the arrival (red line) of the 54.8-80.3 MeV 
protons at the detectors on LI. We also marked the release time 
(solid blue line) as inferred from the VDA analysis and its un¬ 
certainty from the linear fit (dashed blue lines). The dashed blue 
lines mark the “release window”. Eor the 15 January 2005 event, 
we identified type III, II, and IV radio bursts events within the 
release window. 


3. Radio emission associated with SEPs 

3.1. General remarks 

Eor the 83 selected events, we identified the types of radio bursts 
that occurred within each release window. If transient solar radio 
emission (Type III, II, IV-continuum) was present within the pro¬ 
ton release window of an SEP event, we categorized it as an SEP 
with “inferred radio association”. In 18 cases no transient so¬ 
lar radio emission was observed within the release window and 
these events have been characterized as SEPs with “no inferred 
association” (see example in Pigure|2]l. A further examination of 
these events showed that significant radio emission was evident 
relatively near (+30 minutes) the release window in all cases. 
In Pig.|2](f) we present the composite radio spectra of a “no in¬ 
ferred association” case for the 20 August 2002 event. In this 
event there was a strong group of type III bursts 25 minutes after 
the end of the interval ascribed to the release window. We be¬ 
lieve that the inferred lack of association is caused by failures in 
the VDA due to the incorrect determination of onset times from 
the Poisson-CUSUM method. 


2.3. Composite radio spectra 

To compare the proton release time with the associated ra¬ 
dio emission for the 83 selected events directly, we assem¬ 
bled composite dynamic radio spectra. The high-coronal to 
interplanetary (IP) parts of the dynami c spectra consisted of 
Wind/Waves data dBougeret et al.llT995l) with one minute res- 


3.2. Cases with inferred radio association 

After the rejection of the 18 “no inferred radio association” 
events we classified the remaining 65 cases according to the ra¬ 
dio emission within the release window. Type III bursts were 
registered in 53 cases, type II bursts were registered in 44 cases 

* http://secchirh.obspm.fr/ 
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20040725 Wind/Waves * RSTN Dynamic Spectrum 



14:00:00 15:00:00 16:00:00 17:00:00 16:00 



03:30:00 04:00:00 04:30:00 05:00:00 05:30:00 

19980616 Wind/Waves Dynamic Spectrum:18:00:00-21:00:00 



20020417 Wind/Waves * RSTN * ART. Oynam. Spectrum 


. .. 

Mil!* 


IT . - 4 

1 




(e) 


20020820 Wind/Waves ♦ RSTN Dynamic Spectrum 





08:00:00 08:30:00 09:00:00 09:30:00 10:00:00 10:30:00 11:00 


07:30:00 07:45:00 08:00:00 08:15:00 00:30:00 08:45:00 09:00:00 


Fig. 2; Radio spectra for different radio-associated cases: (a) event of 18 June 2000 with type III/II radio association; (b) event of 
11 April 2004 with type III association; (c) event of 25 July 2004 with type III/IV-c association; (d) event of 16 June 1998 with 
type II association, (e) event of 17 April 2002 with type IV-c association; and (f) event of 20 August 2002 with “no inferred” radio 
association. In each panel, the vertical lines have the same meaning as in Fig. [T] 


and type IV-c bursts in 15 cases. In Table IB. 11 we give the proton 
release time and the apparent path length with their uncertain¬ 
ties, as inferred by the VDA for these 65 cases with “inferred 
radio association”; we also give the radio associations for every 
event. 

We further separated the radio associations into more groups 
according to the combination of radio bursts that occurred within 
the release window. For example, if a type III and a type II ra¬ 
dio burst both occurred within the release window we registered 
the event as “type III/II associated” or if only a type III occurred 
within the release window, we registered the event as “type Ill- 
only associated”. In Fig.|2]we present the composite radio spec¬ 
tra for different cases of radio association. For the resulting radio 
associations we specifically found: 25 cases for III/II; 18 cases 
for type Ill-only; seven cases for type Il-only, and 15 cases for 
the mixes of IV-c with III or II. The association rate for type 
Ill-only events is lower than the rate of type III/II group and the 
association rate for type Il-only cases is significantly lower than 
the type Ill-only group. A synopsis of the above radio associa¬ 
tion rates is presented in Table |3 

We also investigated the radio emission associated with the 
65 events without taking the VDA error bars into account. In 
every SEP event we registered the radio emission within +1 hour 
from the proton release time; this range is always larger than the 
VDA errors. With this selection we have a general view of the 
radio burst occurrence, near the proton release time, of SEPs. 
We found that in ten cases we have type Ill-only occurrences; 


Table 2: Overview of the SEPs’ Radio Associations. 


Radio Association Cases: 


III/II 

III 

II 

mix/IV-c 

Total 

No Assoc. 

25 

18 

7 

15 

65 

18 

38% 

28% 

11 % 

23% 

100 % 


t 34 

10 

0 

21 

65 

0 

52% 

16% 

0 % 

32% 

100 % 



Radio association rates within proton release window. 
Radio association rates within +1 h. from proton release time. 


in 34 cases type III and type II; and in 21 cases mixes of type 
IV-c occun-ence (i.e. IV-c, III/IV-c, II/IV-c, III/II/IV-c). Prom the 
association rates, without the use of the VDA errors, we note 
that the type III/II occurrence is still the highest (see Table |2]i. 
Por the cases with only type III occurrence we note that their 
percentage survives, albeit smaller, even if we consider a time 
window of +1 hour from the release time. All the type II events 
are accompanied by type III radio bursts within the time window 
of +1 hour. 

W ith the above analy sis we confirm the well-established re¬ 
sult of ICane et al.l (l2002li that SEPs above ~ 20 MeV are accom¬ 
panied by type III bursts (in our case 53/65 events; i.e. 80%). 
However, from our analysis it is also evident that among these 53 
type III associated events there are 33 cases that are also accom- 
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panied by a type II radio association (33/53 events; i.e. 62%). In 
general, type II bursts were registered in 44/65 cases (i.e. 68%). 

3.3. Flare- or shock-related particle release process? 

From the separation of the radio associations into groups, we 
have concluded that a clear-cut distinction between the flare- 
related and CME-related SEP events is difficult to establisfi. In 
several cases the proton release is accompanied by both a type III 
and a type II radio burst, which are tracers of flare- and shock- 
related particle release processes, respectively. The above con¬ 
clusions apply if we assume that the radio signatures of both 
flare- and shock-related particle release processes inside the re¬ 
lease window can both contribute to the proton acceleration. 

In the cases of type III/II association it is hard to tell if one 
or both related processes can contribute to the particle release 
process. However, we investigated whether the association of the 
proton release times with type II radio bursts indicates a separate 
acceleration process of the SEPs. The results of our investigation 
appear in Section [33] 

3.3.1. Radio spectral characteristics 

Eor the type III/II associated events we further examined the SEP 
composite radio spectra and we determined whether the type III 
bursts emanated from the type II bursts. Eor the cases in which 
this is true, we could argue that the shock wave accelerates both 
the type III emitting electrons and the SEPs. We have considered 
that the type III radio burst emanates from the type II in the cases 
where type II start precedes or coincides with the start of the type 
III. We have not included in the above group the cases where the 
difference between the type III and type II start is less than ~4 
minutes (e.g. the case in Pig.|2la). 

Prom this analysis we have found that in 18 cases the type III 
emanates from the type II (Group 1), while, in seven cases it does 
not (Group 2). Additionally, for the cases with type III/II/IV-c as¬ 
sociation we found that five events fall into the first group, while 
three events comprise the second group. In the majority of the 
events we consider, we did not And the classical sequence of ra¬ 
dio emission in which the type II bursts follow the type III radio 
emission. This result provides a good argument in support of the 
scenario that the type III emitting electrons are accelerated by 
the shock wave. In Pig. [3] we present the composite radio spectra 
for different type III/II radio-associated cases and we show two 
events (Pig.[3a-b) for which the type II radio burst precedes the 
start of the type III and one event (Pig.[3]c) that probably does 
not show this temporal pattern. 

3.3.2. Proton energy spectra 

We have also examined whether the type III/II and Il-only as¬ 
sociated SEP events have a flatter/harder energy spectrum than 
the type Ill-only events. Eor the 65 events with “inferred ra¬ 
dio association”, we have constructed their proton energy spec¬ 
trum by taking the peak flux in each energy channel and pre¬ 
senting it as a function of the proton energy. If we ignore prop¬ 
agation effects then we can argue that the produced peak flux 
spectrum is representative of the injection spectrum (iLin et all 
Il982tl&ucker et a02007h . Although, the peic times occur later 
for lower energies due to velocity dispersion, the peak flux spec¬ 
trum would not be affected from this dispersion. Additionally, 
the calculated spectral slopes from the peak flux spectrum cannot 
depend on the time interval where the calculation is performed. 


We have not included the SOHO/ERNE low-energy channels 
(1 to 10 MeV) in the proton energy spectra analysis to exclude 
events with streaming-limited peak intensities that could cause 
a “plateau” in the energy spectra and also avoid any contamina¬ 
tion that may originate from solar sources that are not related to 
the SEP event under study. In all the cases, the proton energy 
spectra showed a break in the power law with a steeper (softer) 
component at higher energies; this break in the spectrum occurs 
around 25 MeV. The calculation of the spectral power-law in¬ 
dices has been performed by fitting two power-law functions to 
the observed energy spectra, one to the low- and one to the high- 
energy component. This break up in the spectrum can be either 
ascribed to a direct signature of the acceleration process or to 
a secondary process, such as wave-particle interactions, particle 
escape from the acceleration region, or transport effects. 

Eor the derived spectral indices (y), we found values within 
the range -1.27 to -4.99 for the low-energy component (EEC) of 
the spectrum and within the range -1.67 to -4.99 for the high- 
energy component (HEC) of the spectrum. We separated spec¬ 
tral indices into groups according to the radio association of each 
SEP event (see Sect.|3]2). In Table|3]we present a synopsis of the 
mean values for 7 lec and 7 hec for every separate radio associa¬ 
tion group. The separation of y into the above radio association 
groups showed that, in the cases with type II association (i.e. 
III/II, III/II/IV-c, II/IV-c, Il-only), the mean y values are higher 
than the y mean values for the cases with type III association 
without a type II burst (i.e. Ill-only, III/IV-c). 

To validate whether the differences between the derived 
group means are significant, we performed an analysis of vari¬ 
ance (ANOVA) for both the yLEC and yHEC values. The one-way 
ANOVA compares the means between the groups we are inter¬ 
ested in and determines whether those means are significantly 
different from each other. We excluded from the analysis the 
cases of III/IV-c and II/IV-c association because the statistical 
sample was small. Erom the ANOVA we found that the differ¬ 
ence of mean 'Ylec between the type Ill-only and type III/II 
groups is statistically significant (p-value = 2 x 10“^). Similar 
results were found from the comparison of mean ynEC values 
but with higher p-values. Therefore, the spectra of the groups 
with type IIs are harder than those with type Ill-only bursts and 
this result is statistically significant. This is an indication that the 
shock related type II is contributing to the SEP acceleration. 

Moreover, we compared the mean y values between the 
groups of type III/II, Il-only, III/II/IV-c association. We found 
that the ANOVA did not show a clear separation for those groups 
(p-value = 0.18 for EEC and 0.08 for HEC); all the events with 
type II association could belong to the same population of spec¬ 
tral characteristics. The type Il-only group mean y values should 
be closer to a shock related proton acceleration scenario with¬ 
out a flare related counterpart. However, the determined mean 
y values for the type Il-only group was found to have relatively 
high value and this might be explained by the fact that we failed 
to And any type Il-only cases that are completely isolated from 
type III radio bursts, i.e. cases in which the only radio burst in 
the vicinity of the proton release time is a type II burst. We con¬ 
clude that the spectral hardening of the type II cases compared 
to the type III shows that the shock related type II is contributing 
to the SEP acceleration. 

3.3.3. Maximum SEP intensity 

Eurthermore, we examined the maximum SEP intensity in re¬ 
lation to the radio association cases. Earlier studies showed a 
relation between the peak SEP intensities and the associated 
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Fig. 3: Radio spectra for different type III/II radio-associated cases: (a,b) events of 22 November 1998 and 19 December 2002 for 
which the type II radio burst precedes the start of the type III; in these cases the type III emanates from the type II; (c) event of 18 
June 2000, the type III precedes the type II radio burst and we cannot argue that the type III emitting electrons are accelerated by 
the shock wave. In each panel, the vertical blue lines have the same meaning as in Fig[T] 


Table 3: Overview of the SEPs’ energy spectrum indices. 


Radio Association Cases: 



III/II 

III 

II 

III/II/IV-c 

III/IV-C+ 

II/IV-C+ 

Tlec 

-2.41 

-3.22 

-2.79 

-2.35 

-4.25 

-2.28 

Thec 

-2.98 

-3.54 

-3.60 

-2.84 

-3.75 

-2.77 

T t 

^max 

0.043 

0.023 

0.005 

0.111 

0.002 

0.106 


Associations with small statistical sample. 
Maximum SEP intensity in the energy range 54.8-80.3 MeV. 


CME properties, such a s CME speed {e.g. iKahler et al.ll200ll : 
iGonalswamv et alJ2004l) . The correlation between SEP intensity 
and CME speed has been considered strong evidence supporting 
the scenario that SEPs are accelerated by CME-driven shocks. 
We expect that the type II associated SEP events would be more 
intense than the events without a type II association if we assume 
that the type Il-related shock is contributing to the SEP acceler¬ 
ation. However, the scatter of the maximum SEP intensity could 
vary over 3 to 4 orders of magnitude; any correlation with max¬ 
imum SEP intensity would have a high uncertainty. 

For events with “inferred radio association”, we determined 
the maximum proton intensity (Imax) in the energy range from 
54.8 to 80.3 MeV. We found that the events with type III/II asso¬ 
ciation exhibit higher values of I^ax compared to the events with 
type Ill-only association. The group with type III/II/IV-c associ¬ 
ation has the highest value of Imax , meanwhile the type III/IV-c 
showed the lowest value. In Table [3] we present the mean val¬ 
ues of Imax for every radio association group. We further tested 
the statistical signihcance of the Imax mean value differences be¬ 
tween the radio association cases. The ANOVA showed that the 
differences of Imax mean values for all the cases are not statisti¬ 
cally signihcant; therefore, we cannot conclude whether the dif¬ 
ferences originate from the presence (or absence) of a shock re¬ 
lated particle release process; Le. for the type III/II and the type 

III- only group we found a very high p-value {p - 0.45). The 
large scatter of Imax values, which probably originate from parti¬ 
cle transport effects, certainly affect the signihcance of the Imax 
difference between the groups. 

3.4. Acceleration in post-CME current sheet 

We further examined the composite radio spectra for the 15 type 

IV- c associated events (see Sect. |3]l, and investigated whether 
moving type IV bursts were present. For the cases that moving 


type IV bursts could be detected, we could argue that the ener¬ 
getic particles could have been ac celerated at the re connecting 
current sheet behind the CME (see iKlein et ^1201 4ll . 

Erom our analysis we found seven cases (out of the 15 cases 
with type IV-c association) with moving type IV bursts. Eor these 
events, the particle acceleration could have taken place at the 
post-CME current sheet. The most characteristic cases, associ¬ 
ated with moving type IV bursts, are the events of 25 July 2004 
(Eig.l^c) and 17 April 2002 (Eig.l^e). In those events the mov¬ 
ing type IV burst consisted of a prolonged radio emission with 
prominent drift from centimetric to decametric waves. In sum¬ 
mary, we found that in almost half of the cases with type IV-c 
association the particles could have been accelerated at the post- 
CME current sheet. In six out of the seven cases with moving 
type IV bursts, the relative contribution of particle acceleration 
at the post-CME current sheet is difficult to evaluate because the 
moving type IV burst coexists with other transient radio features 
inside the release window. 


4. Locations of the SEP-related flares 


Statis t ical analysis of SEPs (e.g. Cane et al] 19^ ICliver et"^ 
l2004t iKurt et aP l2004t IGonalswamv et al.l l2008h have shown 
that the solar activity associated with SEPs most often occurs 
in the western hemisphere. Erom the distribution s of the “source 
longitu de” of the associated flares lReamesl(ll999ll and iNitta et al.l 
(l2006l) concluded that source regions of the impulsive SEP 
events are distributed about the longitude of best magnetic con¬ 
nection to the observer. The small spread in the longitude dis¬ 
tribution could be attributed to changes in connection longitude 
resulting from variations in solar wind speed. A large scatter in 
longitude away from the nominal P arker spiral con nection has 
been observed for the gradual events. [ReamesI (1 19991) has shown 
that the distribution of the associated flares longitudes for grad- 
ual events is nearl y uniform across the solar disk. Furthermore, 
iKlein et ^ (l2008l) showed that there are cases where open flux 
tubes diverge rapidly with height and these flux tubes connect 
the parent active region to the Parker spiral at the source sur¬ 
face even when the active region is as far as 50° away from the 
nominal connection longitude. 

For the SEP events with “inferred radio association’,’ we in- 
vestigated the loc ations of the asso ciated flares (see T able 3 of 
IVainio et al 120 131 which is based on ICane et al.l (120101) '). In most 
cases (45/65) the associated flares were located in the western 
hemisphere while there was a significant lack of events in the 
eastern hemisphere (5 events only). In 15 cases we did not And 
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Flare longitude (deg.) 

Fig. 4: Distribution of the SEP associated flare longitude in the 
cases with inferred radio association. The red dash-dotted and 
blue dashed lines correspond to different well-connected mag¬ 
netic flux tubes, according to Parker’s model of the interplane¬ 
tary magnetic field, with their longitudes calculated from Eq. |2] 
using different solar wind velocities (i/sw)- The red lines corre¬ 
spond to 50°-60° western longitude (i;sw~415 - 500 km/s). 
The right blue dashed line corresponds to longitude cpo - 80° 
(~312 km/s) and the left blue dashed line corresponds to 
(po - 30° (~ 832 km/s). We also label the 50° -60° longitudes as 
“Well-Connected Region”. 

an associated flare apparently because the events originated from 
far behind the west limb. Both groups with type Ill-only associ¬ 
ation and type III/II had six back-sided cases each; this means 
that there are cases where the radio emission at higher frequen¬ 
cies could be occulted by the limb. 

In Eig. 0] we present the distribution of the associated flare 
longitudes. The maximum of the distribution was found in the 
well-connected region of 50°-60° western longitude, but there 
is a significant scatter on either side of this region from -30° to 
120°. If we accept Parker’s model for evolution of the interplan¬ 
etary magnetic field and a solar wind velocity ~ 450 km/s, 
then the interplanetary magnetic flux tubes have solar connec¬ 
tion longitudes near 50° - 60° western longitude. The equation 
of the solar connection longitude to the conne ction region (i.e. 
Earth ) is given by the following equation (e.g. iNolte & Roelofl 
Il973h : 


fsw r 

Here r is the distance from the Sun centre to the connection point 
(in our case, 212 Rq for the Sun-Earth LI point), a> is the angular 
speed of solar rotation {i.e. the equatorial solar rotation period of 
24.47 days or cu = 2.97 x 10“® rad/sec), Rg - 2.5 R© the heliocen¬ 
tric distance beyond which the magnetic field is assumed to take 
on the form of the Parker spiral, and i/jw is the radial solar wind 
speed, assumed constant in space. In Fig. 01 we label the well- 
connected region of 50° - 60° western longitude with the two red 
lines, calculated from Eq. |2] for solar wind velocity, lijw ~ 415 - 
500 km/s; the blue lines correspond to t/sw ~ 832 km/s (0o ~ 30°) 
and Vsw ~ 312 km/s (0o ~ 80°). 

We further separated the SEPs into groups, according to the 
radio associations inferred from the analysis of Sect. 13.21 We 
separately analysed the cases in which we had type III associa¬ 
tion without type II (Group A); the cases with type II association 
without a type III (Group B) or the mixed cases with type III 
and type II association (Group C). We found that, for the first 


group, the associated flare longitudes were only located in the 
western hemisphere between longitudes 30° to 80° , for 13 out 
of 14 cases. Eor the second group, we found that in four out 
of eight cases the flare occurred within the range -15° to 30°, 
in two cases at longitudes from 30° to 80°, and in two cases 
at longitudes from 80° to 120°. Lastly, the analysis of the third 
group showed a more mixed situation; in nine cases, the flare 
occurred inside the region of 30° to 80° western longitude; in 
another nine cases, at longitudes from 80° to 120°; and, in six 
cases, at longitudes from 0° to 30°. Also, three cases were lo¬ 
cated in the eastern hemisphere. The flare locations of the type 
Ill-related cases are more concentrated in the well-connected re¬ 
gion, whereas the flare locations of the type Il-related cases are 
more scattered from this region. 

5. Proton and electron release timing 

5.1. Proton and type III timing 

Eor the SEPs with inferred radio association, we extended our 
study to the timings between the radio emission of type III and 
the proton release time. Considering that type III radio emission 
is a seconda ry product of Langmui r waves produced by electron 
beams (e.g. iLin et al.lll98lL Il986l) . we used type III bursts as 
proxies of the initial electron acceleration and escape. Prom the 
composite spectra we determined the start time of the associated 
type Ills {i.e. inside the release window). We defined the start 
of a radio burst type III as the point where the radio emission 
was above the background level by three standard deviations, 
preferably at the highest frequency, given that the type III would 
extend to the interplanetary space. In Col. 8 of Table IBTI we 
present the start time of the associated type Ills. 

We found that in all but 12 events, the type Ill-related elec¬ 
tron release occurred before or at the time of proton release 
(time differences from 0 to -24 minutes). The timing distribution 
(Pig. |5^) showed a maximum at -10 minutes with the electron 
release preceding the time of proton release and a mean value at 
-6.4 minutes. We further separated the events according to their 
exact radio associations and examined their timing distributions 
with proton release time. The best synchronization was exhibited 
by events that showed type III/II emissions, with a mean delay 
-5.3 minutes. Eor the events with only type III association we 
found a mean delay -5.6 minutes between type III and proton 
release. The worst synchronization was exhibited by events that 
showed an association with type IV-c radio bursts. 

To validate the results from the proton-type III timing his¬ 
togram, we implemented a separate analysis taking the com¬ 
puted release time un certainty ^t in t o acco unt. We follow the 
method presented by iKrucker et al.l d 19991) to distinguish be¬ 
tween simultaneously released events within the statistical un¬ 
certainties when (f/// - trei)ISt> -1, and delayed events when 
{tin - trei)ISt<-l. Prom the above separation, it follows that 
when the proton release occurred before the type III start, these 
events could only be characterized as simultaneously released; 
in the case where (f/// - trei)l6t > 1 according to the radio associ¬ 
ation analysis (see Section |2l these events have been character¬ 
ized as SEPs with “no inferred association". Eor 31 cases (59%), 
the type Ill-related electron release occurred well before the time 
of proton release, and therefore, these events were delayed; for 
22 cases (41%) the type Ill-related electron and proton release 
was simultaneous. In Pig. [5^, the solid black line denotes the 
smoothed distribution function of the time difference between 
type III start and the proton release time using all the events, 
and the blue and red solid lines denote the distribution func- 
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Type III Timing (All Radio assoc.) 



T<lll> - T<Protons> (min.) 



Fig. 5; (a) The distribution of the time difference between type III 
start and the proton release time derived from the VDA. With the 
arrows we denote the mean and maximum value of the distribu¬ 
tion and with the horizontal lines we mark the errors as inferred 
from the Monte Carlo analysis. The solid lines show the distribu¬ 
tion functions of the type Ill-proton release timings using all the 
events (black) and the separate cases of delayed (red) and simul¬ 
taneous events (blue), (b) The best-fit Gaussians of the type III- 
proton release timings mean and maximum values (black solid 
and dashed lines, respectively) are derived from the Monte Carlo 
analysis using all the events. The red and blue solid lines show 
the best-ht Gaussians to the distributions for the mean values of 
the simultaneous and delayed cases, respectively. 


tion for the separate cases of simultaneous or delayed events. 
The smoothing of the histograms was estimated with a kernel 
smoothing function that returns a probability density estimate 
for a given number of samples. We have scaled the probability 
density estimates for illustration purposes. The timing distribu¬ 
tion of the delayed events (red line) showed a maximum at -9.5 
minutes with the type III start preceding the time of proton re¬ 
lease and a mean value at -10.9 minutes. Additionally, for the 
simultaneous events (blue line), the timing distribution showed 
a maximum at -1 minute and mean value at -0.14 minutes. 

As a consistency check of our timings and to estimate the un¬ 
certainty of the presented mean values, we performed a Monte 
Carlo analysis. With this method we verihed if the above results 
are independent of the selection and the width of the release win¬ 
dow. If we assume that the timing distribution is not sensitive to 
time shifts of the release time, within the error limits of the VDA, 
then if we randomly perturb the release time, the distribution 
characteristics should not change. In our Monte Carlo analysis, 
the perturbation of the measurements was performed by adding 
the release error multiplied by a random number between -1 
and 1 to every proton release time. With this method, we pro¬ 
duced randomly distributed perturbed release times within the 
release error. We repeated the above procedure several times 
(> 10“^ times) to have a statistically significant sample. To anal¬ 
yse the statistical sample produced by the Monte Carlo analysis, 
we tried two different approaches. At hrst, for every “perturbed” 


release sample, we determined the mean value of the time differ¬ 
ences and, from the derived mean values, we analysed their dis¬ 
tribution. Additionally, we applied the above described method 
using the distributions maximum values for every “perturbed” 
release sample. 

In Fig. |5]b we present the results of the Monte Carlo analy¬ 
sis for all the events (black solid and dashed lines for the mean 
and maximum values, respectively) and the separate cases (blue 
and red solid lines for the mean values of the simultaneous and 
delayed cases, respectively). The distributions of the mean val¬ 
ues (solid lines. Fig. |5]b) are localized near the centre line of 
the timing distribution functions and their maximum coincides, 
as expected, with the mean values of the initial timing distri¬ 
butions (Fig. |5]a). This indicates that any randomly distributed 
time shifts of the initial proton release times would produce tim¬ 
ing samples that would have mean values between the limits of 
the Gaussian-fitted functions for each case (black, red, and blue 
htted lines. Fig. |5]b). With the second (maximum) method, we 
found that the distribution is more scattered than that derived 
from the hrst (mean) method and this effect is caused by the 
samples’ discrete form and the histogram binning. In fact, the 
maximum value can signihcantly change between the iterations 
(“perturbations”) and compared with the mean value, which is 
more localized and negligibly affected between the iterations, 
this second approach is less reliable. 

We used the results of the Monte Carlo analysis to match 
the uncertainty of the time differences’ mean and maximum val¬ 
ues. In Fig.|5]a we present the mean and maximum values with 
their error bars for the timing distribution of all the events. The 
size of errors for both mean and maximum in Fig. |5]a are the 
~ +3cr values of the Monte Carlo analysis distributions (con- 
hdence value: 99.73%). For the separate cases, we found that 
the simultaneously released events have a mean time difference 
0.1+5.0 minutes and the delayed events have a mean time differ¬ 
ence -10.9+3.5 minutes. 

5.2. Proton and Wind/3DP electrons timing 

We also analysed the electron data from the Wind/3DP for all 
the type III associated events. In all (53) of these proton events 
we found electron events. We identihed the onset times for seven 
energy channels between 20 keV and 646 keV and we performed 
velocity dispersion analysis using all the available energy chan¬ 
nels. 

For the VDA, we followed the same procedure as described 
in Sect. 12.21 with the exception that we did not use the Poisson- 
CUSUM method for the onset determination. The Wind/3DP 
backgrounds are generally not steady and often show non- 
dispersive variations on various timescales, which can affect any 
statistical method of the onset determination. We determined the 
Wind/3DP onset times by eye, which, compared to the statisti¬ 
cal onset determinations (i.e. Poisson-CUSUM), has the advan¬ 
tage of using the background and rise phase prohle informa¬ 
tion on both short and long timescales (IKahler & Rag(^ l2006t 
iKahler et al .11200^ iHaggertv & Roelof1l2nn2li. 

We performed the VDA and, in most cases (34/53), we 
found a reasonable linear ht. We rejected the events in which 
no dispersion between the Wind/3DP energy channels was ob¬ 
served (15/53) or the path length values were outside the range 
1-3 AU (4/53). A comparison between the electron’s and pro¬ 
ton’s path lengths showed that protons travel considerably longer 
path lengths (s^ = 2.15 AU) compared to the electrons (s^ = 
1.51 AU). The longer path length for protons implies that protons 
are subjected to more scattering from the Sun to the spacecraft 
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Fig. 6: Electron event of 16 October 2000. Top: electron inten¬ 
sity as recorded from Wind/3DP for 7 electron energy channels 
between 20 keV and 646 keV. Bottom: electron velocity disper¬ 
sion analysis. We mark onset times with blue circles and linear 
fit with a red line. 


than the electrons. iMalandraki et alJ (|2012|) also found and con¬ 
firm a similar discrepancy between the path lengths of protons 
and electrons using VDA for a case study. 

In Fig. |6] we present an example for the 16 October 2000 
event. The top panel shows the electron intensity in seven en¬ 
ergy channels. We have added 8.33 minutes to the derived elec¬ 
tron release time to enable the comparison between the arrival 
times of electrons and the electromagnetic radiation at 1 AU. 
For this event, the dispersion between the different energies is 
pronounced. From the electron VDA (Fig. |6] bottom), we calcu¬ 
lated the electron release time at 07:02 UT with an uncertainty 
of +5 minutes, and we found an apparent path length of 1.50 AU 
with an uncertainty of +0.29 AU. 

To compare the proton release from SOHO/ERNE and the 
electron release from the Wind/3DP we plotted the distribution 
of their time difference (Fig. Ha). From the characteristics of 
the distribution, we note that there is a time delay of ~1 min¬ 
utes of the electron release times with respect to the proton re¬ 
lease times. We found that the electron release with the use of 
the Wind/3DP data occurred in most of the cases (33/34, ~91%) 
after or at the time of the proton release. Specifically, we found 
a mean value of 7.1 minutes and a maximum at 5.0 minutes (be¬ 
tween -2.5 and 12.5) for the electron-proton release time dif¬ 
ference distribution (see the solid arrow in Fig. |7]a). We fur¬ 
ther separated the electron-proton release time differences into 
simultaneously released events and delayed events following the 
method presented in Sect lO for the separation of the proton- 
type-III time differences. The uncertainty of the relative timings 
would be, in this case, the sum of the proton/electron release un¬ 
certainties in quadrature, 5t - + dtp, where and is the 

electron and proton release uncertainty, respectively. From the 
separation of the timings into simultaneous and delayed events 
we found that in 18 (~53%) events the proton and electron re- 
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Fig. 7: (a) The grey histogram shows the distribution of the time 
difference between electron release derived from the VDA of 
Wind/3DP electron data and the proton release time, derived 
from the VDA of SOHO/ERNE proton data. The dashed his¬ 
togram shows the distribution of the time difference between 
electron release derived from the VDA of Wind/3DP electron 
data and the type III start. With the solid line arrow and the 
dashed line arrow we mark the mean values of the electron- 
proton and the electron-type III timing difference, respectively. 
The solid lines show the distribution functions of the type III- 
proton release timings using all the events (black) and the sepa¬ 
rate cases of delayed (red) and simultaneous events (blue), (b) 
The best-fit Gaussians of the electron-proton release timings 
mean and maximum values (black solid and dashed lines) de¬ 
rived from the Monte Carlo analysis using all the events. The 
red and blue solid lines show the best-fit Gaussians to the dis¬ 
tributions for the mean values of the simultaneous and delayed 
cases, respectively. 


lease process is simultaneous, while for 16 (~47%) events the 
electron release is delayed compared to the proton release. 

To validate whether the above time differences is indepen¬ 
dent of both proton and electron VDA error we performed Monte 
Carlo analysis in the same way as described in Sect. 15.11 In this 
analysis we only changed the bin size of the maximum distribu¬ 
tion because the statistical sample was small. We present the re¬ 
sults of the Monte Carlo analysis in Fig.|7]b., where it is clear that 
the above timing results are independent of the proton-electron 
timing errors induced by the VDA fit. 

To better illustrate the relative time difference of Wind/3DP 
electron and type III release, we selected to over-plot in Fig.|7]a, 
the relative timing between the electron release from Wind/3DP 
and the start of the type III radio burst, presented with dashed 
line. We note that the electron release occurs on average ~12.3 
minutes after the start of the type III. For the events in which the 
Wind/3DP electron release occurred after the start of the proton- 
associated type Ills (33/34), we found that in all but three cases, 
the durations of the type Ills (mean ~4 minutes) were smaller 
than the electron-type III start time differences. If we addition¬ 
ally consider the uncertainty of the electron release then only 
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six events could be characterized as simultaneously released 
with the type III {i.e. 28/34 were delayed, 82%). Therefore, the 
Wind/3DP electron release typically occurs well after the end of 
the proton-associated type Ills. 

To further investigate other possible implications that can 
be derived by the proton-electron relative timings, we present 
in Fig. [8] the proton-type-III time difference versus the electron- 
proton time difference. With this representation we can separate 
the cases into groups according to the sequence of the type III- 
proton-electron occurrence. For example, the events in which the 
type III occurs first, followed by the proton and the high-energy 
electrons release are located in the labelled region “(Ill/p/e)”. 
Most of the events (19/34) are located in this region, which 
means that type III proton and electron release most probably 
occurs in this successive order. However, there are also some 
cases that do not fit with the Ill/p/e occurrence scenario, for ex¬ 
ample, there are cases in which the proton release is delayed 
compared to the electron release (Ill/e/p) or the proton release 
occurs before the type III (p/III/e). An exampl e where electrons 
are de tected prior to protons has been given bv iMalandraki et ^ 
(120121) . For p/III/e cases (seven events), we note that all the 
events have relative time differences between the type III and 
proton release, which are lower than the proton release uncer¬ 
tainty; the same applies to the single case of the Ill/e/p group 
for the relative time differences between the electron and proton 
release. For the remaining combinations, we found one case in 
p/e/III group. 

Taking the 6t uncertainties into account, we also considered 
the separation of the events into simultaneously released within 
the statistical uncertainty and delayed events as discussed in 
Sect. 15.21 In Fig. 0 we represent the proton-type-III simultane¬ 
ously released events with circles and the delayed events with 
triangles; for the separate cases of the proton-electron release we 
use unfilled markers for the simultaneously released events and 
black-filled markers for the delayed events. We label in Fig. [8] 
the four groups as derived from the above separation. Compar¬ 
ing the number of events in each group we found small differ¬ 
ences; the highest number of events (10) was found in the group 
where the protons and electrons are simultaneously released and 
both are delayed compared to the release of the type III emit¬ 
ting electrons. We note that when we add the electron’s release 
in our analysis and taking the 6t uncertainties into account, the 
separation into groups according to the sequence of the type III- 
proton-electron occurrence becomes ambiguous. 

Summarizing the results of the relative time differences of 
type III bursts start, proton release from SOHO/ERNE and elec¬ 
tron release from Wind/3DP (high-energy electrons) we have 
found that the proton, the high-energy electron release and the 
type III start do not show a dominant release sequence. We found 
a simultaneous proton release compared to the high-energy elec¬ 
tron release (53%); additionally we found that both the proton 
(59%) and the high-energy electron (82%) release were delayed 
compared to the type III start. However, not all the events fit 
with the above scenario, which means that the release sequence 
of protons, the high-energy electrons, and the type III emitting 
electrons is more complex than a single release scenario. 

5.3. Electron release associated radio emission 

Eor the 34 SEP events with an electron event observed by 
Wind/3DP and for which we found a reasonable VDA fit, we 
investigated the associated radio emissions with the electron re¬ 
lease time as inferred by the electron VDA. To identify the types 
of radio bursts that occurred within each electron release win- 



T<Protons>-T<lll> (min.) 


Fig. 8: Scatter plot of the proton-type-III time difference ver¬ 
sus the electron-proton time difference. We use circles for the 
proton-type-III simultaneously released events and triangles for 
the delayed events. For the proton-electron release, we use un¬ 
filled markers for the simultaneously released events, and black¬ 
filled markers for the delayed events. The small black cross rep¬ 
resents the point of the mean value for each time difference. We 
used the red horizontal and vertical lines to separate the scatter 
plot into categories according to the labelled cases {i.e. Ill/p/e, 
Ill/e/p, e/III/p, pAII/e, p/e/III, e/p/III). We label the four result¬ 
ing groups, as derived from the separation, into simultaneously 
released and delayed events; we use rectangles and ellipses only 
for illustration purposes. 


dow, we over-plotted the electron release time and its uncer¬ 
tainty from the least-squares fit at the composite radio spectra 
(see Fig. 0 for an example). We registered the associated radio 
bursts for every event following the same procedure as described 
in Sect. [3] In Fig. |9] we present the composite radio spectra for 
two cases: (1) the event of 18 June 2000 where the electron re¬ 
lease time (green lines) is simultaneous compared to the proton 
release (blue lines) and the type III (black line) burst, and (2) 
the event of 13 July 2004 where the electron release time lags 
compared to the proton release and the type III burst. 

From the inspection of the composite spectra, we registered 
radio bursts associated with the electron release for all the cases 
except from the event of 07 October 1997. For this event a pos¬ 
sible reason for the lack of association might be that the relevant 
radio emission was occulted as the associated flare was located 
well behind the solar limb. In the remaining 33 cases with in¬ 
ferred radio association we registered 23 cases with type III and 
II association (III/II); seven cases with only a type III radio burst, 
and three cases with mixed type IV-c association. 


6. Estimate of proton release heights 

For the 44 SEP events associated with type II bursts, we es¬ 
timated the proton release heights using the assumption that 
type Il-related shock wave actually accelerates the first particles. 
Since we have shown in Sect. 13.31 the accelerative role of the 
type Il-related shock waves, the above assumption could not be 
rejected entirely although we cannot determine the exact relative 
contribution of the shock-related and flare-related processes to 
the SEP acceleration. Eor each event, we determined the corre¬ 
sponding CME observed by the SO HO’s Larne Angle and S pec- 
trometric Coronagraph (LASCO: lBrueckner_e^^ 1995) and 
registered on the online catalogue of lYashiro et al.l ( 2004 ). We 
found associated CMEs for all the events, except one case dur¬ 
ing a SOHO/LASCO data gap. In all cases there was only one 
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Fig. 9: Composite radio spectra for two different electron events. 
Top; event of 18 June 2000. Bottom; event of 13 July 2004. With 
the vertical blue line and dashed lines we present the proton re¬ 
lease time derived from the VDA analysis and its uncertainty, 
respectively; with the vertical green line and dashed lines we 
present the electron release time derived from the VDA analysis 
and its uncertainty, respectively; with the black dash-dotted line, 
we present the start of the type III burst. 



Fig. 10; The distribution of the estimated proton release heights 
from the CME leading edge height-time measurements. The blue 
solid line represents the distribution function of the estimated 
proton release heights for the SEP events whose source region 
longitude was greater than 80°. 

CME near the onset of the proton event and specihcally near the 
proton release window. 

Eor the estimate of the proton release heights, we used the 
information for the height of the CME leading edge as mea¬ 
sured by the automated procedure described by lYashiro et alJ 
(I 2 OO 4 I 1 . We extrapolated (or interpolated) their projected height¬ 
time measurements to the time of the proton release as estimated 
from the VDA. Erom the resulting distance, we have a value for 
the proton release height if we accept that the protons are primar¬ 
ily released at the leading edge of the CME. Recent studies have 
shown the importance of the lateral CME expan sion as a driver 
of the particles into well-connected regions (iRouillard et al.l 


2 OI 2 I). Another important issue is the standoff distance (e.g. 
Ontiveros & Vourlidasll200^ iGonalswamv & Yashiroll20Tl1l be¬ 

tween the CME leading edge and the shock front that might in¬ 
troduce an error to our estimates. Actually, the particles are not 
necessarily accelerated at the CME leading edge, but they can be 
located anywhere in the sheath region between the CME front 
and the shock front. Thus, our method might underestimate the 
“real” proton release height, depending on the standoff distance 
and the particles spatial distribution into the sheath region. 

In Eig. [To] we present the distribution of the resulting SEP 
projected release heights. The height is measured from the solar 
centre, so the photosphere is at a height equal to one. In eight 
cases, we found SEP release heights that correspond to heights 
below the photosphere (< 1 R 0 ) and these cases are artefacts of 
the analysis. This is probably an effect of the initial assumption 
that the CME leading edge propagation can be approximated as a 
linear ht of height or result from projection effects since the tra¬ 
jectory of a CME is determined in projection on the plane of the 
sky. Our assumption that the height increases linearly with time 
and usually agrees with CME height-time data above 1.5 Rq, but 
the lack of measurements for heights below the LASCO/C2 held 
of view can, in some cases, affect the accuracy of the computa¬ 
tion. 

The SEP projected release height distribution (Pig.fTOli peaks 
at 3.5 Rq (between 3 - 4 Rq ) and extends up to 8 Rq ; we found a 
mean value of 3.4 Rq, excluding from the calculation the hve 
cases with heights below the photosphere. Eor the cases with 
type Il-only association, we found for the proton projected re¬ 
lease heights a mean value of ~4.4Rq. Eor these cases, we do 
not present the peak of the distribution because the number of 
cases is too small to support a statistically signihcant maximum. 
Moreover, we over-plot in Eig. (TOj the distribution function of 
the estimated SEP projected release heights only for the SEP 
events whose source region longitude was greater than 80°. We 
used this criterion to reduce the sample of events with uncer¬ 
tain height-time prohles due to projection effects. Prom the re¬ 
sulting distribution we found a maximum at ~2.25Ro and a 
mean value at ~3.79Ro. Additionally, the cases with unrealis¬ 
tic heights (<lRo) were signihcantly reduced. Our hndings can 
be interpreted by the following scenario. As a CME expands ra¬ 
dially outward from the solar corona, a shock wave forms when 
the Alfven speed ahead of the CME falls below the CME speed. 
So it is not surprising that the SEP-associated shocks accelerate 
particles at starting heights near 2.5 Rq if we accept the A lfven 
speed-height dependence. According to iMann et al.l (l2003h . the 
Alfven speed decreases rapidly with radial distance to a mini¬ 
mum at ~1.2- 1.8 Rq, then rises to a maximum at ~3.8Rq and 
decreases thereafter. 

Similar results of rel ease heights have been re ported 
in previous st u dies (i Hutt unen- He i kinmaa et al.l l2005[ 
IGonalswamv et all 120121; fReamesn2009h Tor pr oton/helium 
SEPs or ground-level events (GLEsl. lReamesI (l2009l) has studied 
the onset times in large SEPs, which have been detected by 
neutron monitors at ground level, and converted this release 
time to a radial distance of the source shock wave from the 
Sun. They found that the acceleration for well-connected 
ev ents begins at 2-4Rq. A lmost the same result was found 
by IGonalswamv et all (1201 2h for the release of several GLEs; 
they showed that the release occurs when the CMEs reach 
an average height of ~3 . 09Rq (between 1.71 to 4.01 Rq). 
IHuttunen-Heikinmaa et al.l (l2005h also determined the apparent 
proton release heights for 13 proton/helium events with the 
most reliable release times (path length; 1.0-1.5 AU) and they 
found a release height from 2 to 10 Rq. The above proton release 
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heights have been converted to distance from Sun’s centre, when 
necessary, to be comparable with our measurements. 

7. Summary of results 

We have studied the properties of major SEP events, especially 
focussing on (1) the association of the SEP release time, as in¬ 
ferred by the velocity dispersion analysis (VDA), with transient 
solar radio emissions recorded by space and/or ground based ra¬ 
dio spectrographs; and (2) the time difference of proton release 
with respect to the escape of keV electrons into space, derived 
from both DH type III bursts and VDA. Our key findings are: 

- Both flare- and shock-related particle release processes are 
observed in high-energy proton events at >50 MeV. A clear- 
cut distinction between flare-related and CME-related SEP 
events is difficult to establish. 

- Proton release is most frequently accompanied by both type 
III and II radio bursts (38%), but there is a significant per¬ 
centage of cases with only type III occurrence (28%). 

- Typically, the protons are released after the start of the asso¬ 
ciated type III bursts and simultaneously or before the release 
of energetic electrons. 

- The locations of the major SEP-related flares was concen¬ 
trated in the western hemisphere with a peak in the well- 
connected region of 50° - 60° western longitudes and signif¬ 
icant scatter in the range from -30° to 120°. 

- The proton release for the type II associated cases typically 
occur at heights from 2.0 to 3.5 Rq. 

8. Discussion 

Unlike previous statistical studies, we managed to directly asso¬ 
ciate the proton release time, as inferred by the proton VDA, with 
different types of transient radio emission, from the composite 
radio spectra, within the proton release window. Based on our 
radio association analysis, we fail to find support for the hypoth¬ 
esis that only the flare-related particle release process dominates 
in SEP events with energy above ~50MeV. If the latter hypoth¬ 
esis were true we would expect that the proton release should 
mainly be associated with the occurrence of type Ills only. How¬ 
ever, a significant number of cases shows radio association char¬ 
acteristics, which resemble those of the classical group of im¬ 
pulsive cases where the SEPs observed in high energies are only 
associated with flare-related phenomena, such as the presence of 
type III radio bursts. Additionally, we found that, less frequently, 
SEP events are only associated with shock-related phenomena, 
such as the presence of type II radio bursts. Eor the type IV-c ra¬ 
dio associated events, we found that an acceleration of particles 
in the post-CME current sheet could apply in almost half of the 
cases, although the relative contributions from various sources 
is difficult to evaluate because the moving type IV bursts coex¬ 
ist with other transient radio features. In summary, the results of 
the radio association analysis leads us to conclude that, in many 
cases, it is possible that the flare-related and shock hypothesis of 
particle release could both apply. 

The above conclusions apply if we assume that all radio- 
related processes (flare or shock related) inside the release win¬ 
dow can both contribute to the proton acceleration and release. 
In general it is hard to derive the relative contributions of these 
processes to the particle release process, but we investigated 
whether the association of the proton-release time with type 
II radio bursts indicates a separate acceleration process of the 
SEPs. Prom the inspection of the type III/II radio spectra, we 


found that in most of the cases (73%) the type III bursts emanate 
from the type II bursts. This result supports the scenario that the 
type III emitting electrons are accelerated by the shock wave for 
those cases. As an additional argument, we computed the proton 
energy spectra for all the SEP events and found that the events 
with type II association have a harder energy spectrum than those 
without a type II. Summarizing the above results, we have found 
that the presence of a type II may indicate a contributing accel¬ 
eration process from the associated shock wave. 

In Sect. |5] we presented the relative timing between the pro¬ 
ton release and the start of the associated type III radio emis¬ 
sion. We used type III bursts as proxies of the initial electron 
acceleration and escape. In general, we found a late release for 
protons. In 59% of our cases, the type Ill-related electron release 
occurred well before the time of proton release; these cases have 
been characterized as delayed. The apparent delay between the 
proton release and the type III radio emission can be ascribed to 
several reasons, such as selective acceleration, transport effects, 
particle diffu sion, and particle trappi ng. An other possibility was 
suggested bv lKlein & PosneJ (l2005h . who attributed the delay of 
protons to the successive injection of particles in different flux 
tubes, which are not all connected to the spacecraft. Lastly, we 
did not And a dependence between the timings and the flare lon¬ 
gitude, nevertheless, some events outside the well-connected re¬ 
gion exhibit the worst synchronization. 

We further investigated the relative timing between the pro¬ 
ton and electron release using Wind/3DP data. When we im¬ 
plemented the dt uncertainties in the analysis and we sepa¬ 
rated the events into simultaneously released and delayed events, 
we found that in 53% of our cases the electron release (ob¬ 
served by Wind/3DP) compared to the proton release (observed 
by SOHO/ERNE) was simultaneous within the statistical un¬ 
certainty. As far as the remaining delayed cases ar e concerned 
we no te that similar results have been presented bv ICliver et al.l 
(Il982l) . Of course, the com parison of our relative timings with 
those of ICliver et al.l (Il982h might be indirect mainly because of 
the different electron-proton energies analysed in each study, but 
the underlying physical reason could be similar. This delay could 
be attributed to turbulence or waves associated with the CME- 
driven shock that traps the low-rigidity electrons much more ef¬ 
fectively than the protons, selective acceleration, or transport ef¬ 
fects. 

As far as the delayed electron release is concerned, com¬ 
pared with the start of the type III (~12 minutes) we additionally 
found that in all except six cases the durations of the type Ills 
(mean: 4 minutes) were smaller than the electron-Type III start 
time differences. Therefore, the Wind/3DP electron release typi¬ 
cally occurs well after the end of the proton-associated type Ills. 
Previous s tudies have investiga ted this delay providing different 
scenarios. iKrucker et alJ (1 19991) attributed this delay to two dis¬ 
tinct populations of electrons, the low-energy population associ¬ 
ated with the type III radio burst and the high-energy population 
representing the delayed electrons. Some of our analysed cases 
(17/33) could fit into the above scenario because the inspection 
of the composite spectra showed a very good association of the 
electron release with group of type Ills or isolated type III sepa¬ 
rate from the initial type III, which occurs at the start of the SEP 
event. _ _ 

iHaggertv & Roelol (l2002l) showed that the electron delay 
with respect to the type Ills could be due to the acceleration and 
release of the near-relativistic electrons by an outgoing coronal 
shock. However, we did not And any dependence of the re lative 
timing s with the existence (or not) of a type II association. ICanel 
(l2003l) showed that the low- and high-energy electrons belong 
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to the same population and, thus, interaction effects in the inter¬ 
planetary medium might cause the delays. Even though the elec¬ 
tron propagation and the interaction effects in the interplanetary 
medium could be a possible cause for the observed electron de¬ 
lays, we did not find any correlation between the relative timings 
and the electron apparent path lengths. 

In Sect. |4] we examined the locations of the related flares 
for the SEPs with “inferred radio association”. The maximum 
of the flare helio-longitude distribution was found in the well- 
connected region of 50° to 60° western longitude, but there is 
also significant scatter on either side of this region from -30° 
to 120°. We further separated the SEPs into groups, according 
to the radio associations inferred from the analysis of Sect. 13.21 
Prom this separation, it is evident that the type Ill-only cases are 
associated with flare locations that are more concentrated to the 
well-connected region, whereas for the cases where the proton 
release is associated with type IIs only (shock-related cases), the 
SEP-related flares could be located far from the well-conne cted 
region. A similar result has been presented bv lReame^ (Il999h for 
the cases of the SEP associated flare locations for the impulsive 
and gradual cases. 

Por the SEP events related to type II bursts and CMEs (44 
events), we estimated the proton release heights from the height¬ 
time profile of the CME leading edge at the moment of the pro¬ 
ton release as inferred from the VDA. The height distribution of 
proton release peaks at 2.5 Rq (between 2-3 R©) and extends up 
to SRq. Previous studies have prese nted similar results for the 
release heights of SEPs and GLEs d Huttune n-Heikinrna a et af 
l2005t l^palswamv et alj|2008l : [Reamesl]2009ll . Unlike [R eamesi 
(!2009h . we found no dependence of the SEP release height with 
the associated flare longitude. 
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Appendix A: Calculation of energy density of 
energetic protons 

Here we investigate the contribution of the energetic particles to 
the total energy content of the plasma by calculating the energy 
density of a major SEP event at 1 AU and by comparing it with 
the energy density of the local magnetic field. In the case where 
the computed energy density of the energetic protons (ep) at 1 
AU exceeds the local magnetic field energy density (ee), the par¬ 
ticles are no longer confined into magnetic flux tubes and they 
can modulate the local magnetic field. The proton energy density 
can be derived from the following equations: 


r r^^dEdQ. r r^^dlnEdQ, 

Jn Jei v Jn Je, Pc dE 


(A.l) 


where E is the particle energy, Q is the instrument solid angle, 
and p is the particle velocity in it/c. The integration over the an¬ 
gle Q is calculated for the particle detector effective angle (field 
of view) and the integration over the energy is calculated from 
the minimum (£1) to the maximum (£2) energy of the detected 
protons. The magnetic energy density is given by the following 
equation: 


25 



Mev/n 


= [Jhn\ (A.2) 

2£o 

In the above equation, the magnitude of the magnetic field 
is measured in Tesla and the magnetic energy density is in 
Joules/m^. 

We selected to analyse the SEP event of 26 December 2001. 
To derive proton spectra, we used the SOHO/ERNE particle data 
of one minute average values. T he SOHO/ERNEs’ particle de¬ 
tectors field of view (half angle; iTorsti et Slll995h is ~30° for 
the low-energy detector (LED) and ~60° for the high-energy de¬ 
tector (HED). Combining the data of the HED and the LED, we 
have an energy coverage from 1.68 to 108 MeV. We first esti¬ 
mated the SEP energy spectra and then we integrated £ dA/d£ 
to derive the Cp. Eor the calculation of the energy spectra, we 
chose not to perform a time shift between the different energy 
channels because we wanted to compare Cp with the measure¬ 
ments of the local magnetic field, at 1 AU and not at the accel¬ 
eration site. Eor the computation of the cb, we used one minute 
average data from Wind Magnetic Field Investigation (MEI) in¬ 
strument, for the magnitude of the average magnetic field vector. 

In Eig. lA.ll we present the proton intensity and the magni¬ 
tude of the magnetic field for the event of 26 December 2001. At 
12:10 UT we calculated the proton energy spectra (blue line) and 
the measure of £^ dA/d£ (green line). The latter in logarithmic 
scale is proportional to the integral in Ea. IA.ll and its value shows 
the energy region where the contribution to the integral is maxi¬ 
mum; in our case the maximum contribution to the proton energy 
density comes from protons with energy close to 10 MeV. The 
magnetic energy density was eB=2.6 10“^'J/m^ and this value 
is almost three orders of magnitude larger than the calculated 
proton energy density whose value was ep=3.8 10“'^ J/m^. 

Erom the above analysis, it is evident that the magnetic en¬ 
ergy density for the 26 December 2001 event is much larger than 
the non-thermal proton energy density. We repeated the same 
analysis for the very energetic events (Ip > 0.2 pfu) of 26 Octo¬ 
ber and 02 November of 2003 and we found similar results; the 
lowest difference between eB and Cp was almost two orders of 
magnitude. 


Eig. A. 1: Data for the calculation of the proton and magnetic en¬ 
ergy density for the SEP event of 26 December 2001 at 12:10 UT. 
Top: average of the magnetic field vector from Wind/IMP. Mid¬ 
dle: proton intensity in 20 energy channels of SOHO/ERNE. The 
vertical black dash-dotted line marks the time of the energy den¬ 
sity calculation. Bottom: proton energy spectrum in logarithmic 
scale at 12:10UT shown with a blue line and the measure of 
E^dNIdE shown with a green line. 

Appendix B: Table of events with inferred radio 
association 

In Table IB. ll we present our results for the 65 proton events with 
“inferred radio association”. In columns 2 and 3 we give the pro¬ 
ton release times and path lengths, respectively, with their uncer¬ 
tainties as inferred from the VDA. Also, in column 4 we give for 
every event the resulting transient radio emission associated with 
the proton release according to the analysis of Sect. [3] In Col. 5 
we present the locations of their associated flares that were used 
in the analysis of Sect. |4] The entries correspond to flares 
occurred far behind the west limb. The electron release times 
and path lengths with their uncertainties as inferred by the VDA 
of Sect. 15.21 are given in columns 6 and 7, respectively. The esti¬ 
mated start times of the type III bursts used in the relative timings 
of Sect. 15. II are presented in Col. 8. 
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Table B.l; Events With “Inferred Radio Association”. 


Event Day 

Proton 

Rel. (UT) 

Proton Path 
Length (AU) 

Proton Rel. 
Radio Assoc. 

Flare 

Long. 

Electron 
Rel. (UT) 

Electron Path 
Length (AU) 

Type-III 
Start (UT) 

19971007 

12:53 + 07 

2.66 + 0.13 

Ill 

_ 

13:13 + 03 

1.48+0.18 

12:51 

19971104 

06:06 + 06 

1.86 + 0.15 

III/II/[IV-c] 

33W 

06:13 + 03 

1.17 + 0.15 

05:56 

19971113 

21:30 + 06 

2.49 + 0.12 

III 

- 

21:36 + 02 

1.02 + 0.10 

21:20 

19971114 

13:16+10 

2.42 + 0.27 

III 

- 

13:16 + 04 

1.54 + 0.24 

13:00 

19980420 

10:15 + 05 

2.67 + 0.12 

III/II 

90W 

10:25 + 01 

1.15+0.07 

10:03 

19980502 

13:45 + 02 

1.34 + 0.06 

III/II/[IV-c] 

15W 

- 

- 

13:34 

19980509 

03:34 + 06 

2.55 + 0.12 

III/II/[IV-c] 

lOOW 

- 

- 

03:22 

19980616 

18:54 + 07 

1.80 + 0.11 

II 

115W 



- 

19981122 

06:40 + 03 

1.71+0.08 

III/II 

82W 

06:51+06 

1.45+0.28 

06:41 

19990424 

13:21+06 

2.14 + 0.12 

II 

- 

- 

- 

- 

19990509 

18:08 + 05 

1.73 + 0.13 

III 

95W 

18:13 + 03 

1.03+0.15 

17:59 

19990527 

10:41+03 

1.77 + 0.08 

III 

- 

- 

- 

10:38 

19990601 

18:56 + 06 

2.71+0.15 

III/II 

_ 

19:04 + 03 

1.85+0.16 

18:44 

19990611 

00:41+03 

1.69 + 0.08 

III/II 

_ 

- 

- 

00:41 

20000212 

04:28+13 

2.17 + 0.42 

II 

24W 

- 

- 

_ 

20000217 

20:50 + 05 

1.54 + 0.08 

III/II 

07E 

- 

- 

20:29 

20000218 

09:32 + 05 

1.47 + 0.14 

III/II 

- 

- 

- 

09:23 

20000303 

02:13 + 06 

2.03 + 0.16 

III/II 

60W 

02:23 + 02 

1.02 + 0.10 

02:12 

20000404 

15:28 + 03 

1.51+0.07 

III 

66 W 

- 

- 

15:17 

20000618 

01:56 + 05 

1.77 + 0.24 

III/II 

85W 

02:00 + 05 

1.49 + 0.27 

01:57 

20000722 

11:26 + 04 

1.84 + 0.09 

II 

56W 

- 

- 

- 

20000912 

12:07 + 05 

2.46 + 0.14 

III/II/[IV-c] 

09W 

12:21+04 

1.38 + 0.21 

11:46 

20001016 

07:02 + 05 

1.63 + 0.10 

III/II 

95W 

07:02 + 05 

1.50 + 0.29 

06:52 

20001025 

10:11+06 

2.59 + 0.13 

III/II 

120 W 

10:11+05 

2.69 + 0.26 

10:01 

20001124 

05:04 + 05 

2.35 + 0.11 

III/II 

03W 

05:07 + 03 

2.15+0.20 

04:58 

20010105 

17:39 + 06 

2.39 + 0.12 

III 

- 

17:44 + 04 

2.24 + 0.21 

17:33 

20010329 

09:56 + 07 

2.97 + 0.15 

III 

12W 

10:02 + 08 

1.54 + 0.49 

09:59 

20010402 

11:28+16 

2.59 + 0.53 

II/[IV-c] 

62W 


- 

- 

20010410 

05:34 + 05 

2.31+0.11 

III 

09W 



05:13 

20010415 

13:47+14 

1.61+0.16 

III/II 

84W 

- 

- 

13:39 

20010507 

12:37 + 03 

1.62 + 0.08 

II 

_ 

- 

- 

- 

20010520 

06:02 + 06 

2.31+0.13 

III/II 

91W 

06:18 + 01 

1.50 + 0.86 

06:02 

20010604 

16:29 + 05 

1.72 + 0.09 

III/II 

60W 

16:30 + 03 

1.07 + 0.19 

16:22 

20010615 

15:40 + 05 

1.74 + 0.14 

III/II 

- 

15:41+02 

1.49 + 0.13 

15:34 

20010619 

03:42 + 07 

2 .21+0.20 

III/II 

- 

- 

- 

03:23 

20010915 

11:48 + 04 

1.51+0.08 

III/II 

53W 

11:46 + 03 

1.02 + 0.13 

11:48 

20010924 

10:17 + 07 

2.63 + 0.16 

III 

23E 

- 

- 

10:18 

20011019 

01:20 + 09 

2.35 + 0.32 

II/[IV-c] 

18W 

- 


- 

20011019 

16:31+08 

2.55 + 0.20 

III 

29W 

16:49 + 03 

1.75+0.14 

16:30 

20011022 

15:14 + 05 

1.88 + 0.09 

III/II 

18E 

- 

- 

15:02 

20011226 

05:27 + 03 

1.37 + 0.07 

III/II/[IV-c] 

54W 

_ 

_ 

05:13 

20020127 

13:03 + 06 

2.18 + 0.17 

II 

_ 

_ 

_ 

- 

20020220 

05:54 + 08 

1.36 + 0.11 

III 

72W 

- 

_ 

05:55 

20020318 

03:18 + 14 

2.98 + 0.44 

II/[IV-c] 

lOOW 

- 

- 

- 

20020417 

09:27 + 06 

1.84 + 0.12 

[IV-c] 

34W 

- 

- 

- 

20020421 

01:15 + 06 

1.24 + 0.11 

III/II 

84W 

01:20 + 03 

1.37 + 0.17 

01:20 

20020707 

11:30 + 06 

1.87 + 0.10 

III 

95W 

11:39 + 01 

1.00 + 0.07 

11:17 

20020814 

01:57 + 04 

1.36 + 0.08 

III/II 

54W 

- 

- 

02:01 

20020818 

21:38 + 07 

1.46 + 0.25 

II 

19W 

- 

- 

- 

20020822 

01:46 + 06 

2.03 + 0.16 

III 

62W 

02:05 + 06 

1.88+0.31 

01:50 

20021109 

13:13 + 14 

2.99 + 0.32 

III/II 

29W 

13:25 + 01 

2.16 + 0.07 

13:08 

20021219 

21:40 + 04 

1.74 + 0.10 

III/II 

09W 

21:51+01 

1.02 + 0.06 

21:41 

20030531 

02:17 + 05 

1.85 + 0.13 

III 

65W 

02:20 + 03 

1.94 + 0.15 

02:21 
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Table B.l: Continued. 


Event Day 

Proton 
Rel. (UT) 

Proton Path 
Length (AU) 

Proton Rel. 
Radio Assoc. 

Flare 

Long. 

Electron 
Rel. (UT) 

Electron Path 
Length (AU) 

Type-III 
Start (UT) 

20031028 

10:21+07 

2.21+0.16 

II/[IV-c] 

08E 


_ 

— 

20031102 

09:27 + 05 

1.91+0.16 

III/II 

- 

- 

- 

08:56 

20040411 

04:14 + 06 

2.14 + 0.11 

III 

46W 

04:15 + 03 

1.31+0.14 

04:08 

20040713 

00:12 + 09 

2.77 + 0.21 

III/II 

59W 

00:28 + 03 

1.28 + 0.12 

00:14 

20040725 

14:55 + 09 

2.60 + 0.20 

III/[IV-c] 

33W 

- 

- 

14:45 

20041101 

05:34 + 04 

2.00 + 0.10 

III 

_ 

05:50 + 03 

1.45 + 0.15 

05:42 

20050115 

06:23 + 03 

1.67 + 0.08 

III/II/[IV-c] 

06E 

06:23 + 01 

1.38 + 0.09 

06:07 

20050713 

14:15 + 08 

2.77 + 0.18 

III/[IV-c] 

80W 

14:16 + 02 

1.27 + 0.12 

14:02 

20050822 

01:19 + 05 

1.89 + 0.10 

III/[IV-c] 

48W 

- 

- 

01:20 

20060706 

08:27 + 05 

2.71+0.13 

III/II 

32W 

08:37 + 04 

1.55 + 0.21 

08:21 

20061213 

02:21 + 11 

2.13 + 0.41 

III/II/[IV-c] 

23W 

- 

- 

02:24 

20100814 

10:05 + 04 

1.50 + 0.08 

III/II/[IV-c] 

54W 

10:06 + 03 

1.02 + 0.15 

10:03 
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